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INTERACTION OP THE INTERMEDIATE ENERGY NEUTRINO 

WITH NUCLEI 

E. V. Bugayev, M. A. Rudzskly, G. S. Bisnovatyy- 
Kogan, Z. P. Seidov 


1. Introduction 




Many authors have made theoretical studies of the Interaction 
of neutrinos with small and average energies. The research object 
Is usually the transition cross sections In definite states of the 
finite nuclei, and also the total (with respect to all of the nu- 
clear excitations which are possible for a given neutrino energy) 
cross sections of the v -capture. It Is usually difficult to cal- 
culate the cross section In the energy region E^'*'I0*^I00 MeV, 
where a knowledge of the structure of the low-lying excited states 

of the target nucleus Is very Important. Therefore, at present 

12 

detailed calculations only exist for the simplest C [1], [2], 

^^0 [3], [^], ^"^Cl 151 , [ 6 ] nuclei. In the energy region which Is 
far from the threshold (E^ > 100-200 MeV), it Is possible to use 
simple universal models, such as the model of the Permi-gas [7] 
and the envelope model of J. Bell and S. Llewellyn Smith (BLS) 

[8J, which use the approximation of completeness. 

Further research on the neutrino-nuclear reactions Is impor- 
tant for two reasons. 

1 , The low energy neutrinos play an important role In astro- 
physics. In the case of steady state combustion, particularly In 
the case of a collapse, the stars emit neutrinos which can be re- 
corded on the Earth. The energy of the emitted neutrinos In the 
case of steady state combustion for stars of the main sequence 
does not exceed 8 MeV, and It may be much greater in the case of 
a collapse. In the first estimates made by Ya. B. Zel'dovlch and 
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0. Kh. Guseynov [ 9 ], the energy of neutrinos emitted in the case 
of a collapse was 60 MeV. In the recent calculations of D. K. 
Nadjozhin, who made the most detailed examination of neutrino 
processes in the case of collapse, it was found that in the case 
of diffusion in the nucleus of a star which is not transparent 
for neutrinos their energy decreases and Evinax ~ 20-30 MeV [10] . 

The cross section of the V -capture must be known not only to 
study the possibilities of recording the neutrinos, but also tb. 
study the collapse dynamics. We must also stress that, in princi- 
ple, there may be fluxes of much higher energy neutrinos having 
an astrophyslcal origin, for example, neutrinos from pions which 
decompose in a year, etc. Thus, the neutrino energy region from 
the threshold up to approximately 100 MeV is of interest for 
astrophysics. In this region, as was already noted above, the 
calculations are very complex and are made for each nucleus in- 
dividually, Simple nuclear models close to the threshold are 
not accurate enough, especially for light nuclei. 

2. It is known that nuclear targets may be Widely used to 
Increase the statistics in accelerator experiments . If the aver- 
age neutrino energy is not very large (for example, as for the 
Argonne accelerator, where ^ 200-300 MeV), then a detailed 

study is needed of the effects of the nuclear structure. Within 
the framework of the shell model of the nucleus, this problem has 
been studied in [8] . The authors used wave functions of a harmonic 
oscillator for the calculations, so that they could consider the 
continuous spectra of the nucleus excitations. The excitation 
spectrum in this model is far from the real one, and it is diffi- 
cult to establish the error connected with it. The most detailed 
calculations at present have been made for the v-capture react.i ons 
in and in [1-4] . They obtained the total cross sections 

in the neutrino energy range from zero to 500-600 MeV. The article 
lU studied the v-capture by the :nucleus Experimental data 

were used from 22 low-lying excited states in obtained from 

experiments on the electro- and photoexcitation of Similar 

calculations were made for the nucleus 131 . Calculations 


2 


were made in [2] of the cross section of seml-lepton processes in 
nuclei with A = 12 on the basis of the shell model of the nucleus 
with an oscillator potential hole and allowance for the residual 
interaction. For agreement with the experimental data, all the 
transition cross sections had to be reduced. In spite of consider- 
ing the residual interaction. A similar calculation was made for 
in [4] . 

Detailed calculations were made In [1-4] and are only possi- 
ble for very light, well-studied nuclei. It is extremely complex 
to perform calculations for heavier nuclei, for example, for ^*^01 
[6], since the number of low-lying excited states in the nuclei 
rapidly Increases with an increase in their atomic weight A. 

The purpose of this article is to formulate a unified method 
for calculating the cross sections of inelastic neutron=nuclear 
reactions for the neutrino energy regions from threshold values to 
approximately 300 MeV. A unified description is given in this 
region using the precise calculation of all the nuclear transitions 
described within the framework of a single-frequency shell model 
with the Wood-Saxon potential CWS) . The WS potential has a dis- 
crete and continuous spectrum, and gives the correct sequence of 
single particle levels up to the very heavy nuclei. Since present- 
ly it is impossible to give a precise description of the true 
nuclear states within the framework of this model, to obtain an 
approximate total transition cross section in all the possible 
finite states of the nucleus, vie calculate the total transition 
cross section in all the finite single particle states, in which 
the transitions obey the law of conservation of energy. In real 
calculations, certain simplifications are made, noted in the 
second section. 

We consider elastic and inelastic reactions with the neutrino 
electrons : 


( 1 . 1 ) 


♦ (A.l) e- * 

Preliminary results of the calculations were published in [11]. 

This article consists of four parts and the appendices. The 
second part gives the derivation of the formulas for the rross sec- 
tions of the reactions (1.1) - (1.3). The third part discusses the 
results of calculating the cross sections and makes a comparison 
with the work of other authors. In conclusion, the methods of im- 
proving this procedure are discussed. Certain computational details 
are examined in the appendices . 

2 . Interaction cross section in single particle shell model 


1. Charged flux 


The Interaction cross section for a neutrino with a nucleus, 
ly transitions in an arbi’ 

•s ' (in a laboratory system): 

\ "n -1* ) ci-k. 

( 2 . 1 ) 


accompanied by transitions in an arbitrary group of nuclear states 
j. equals 

^ ^ ' oL k' 






Here, k,k,' and k*. 4 - s.re the pulses of the neutrino and electron, 
respectively; q = k-k , and — energy of the initial and 
final nuclei state; — lepton tensor, 

1 , ± kjk', ) ^ (2.2) 

Here 3. is the 4-vector, determined by the expression 

r 

*^We use the system of units • C • i , and the metrics is such 
that 


where lt,> and are the final and Initial nuclei states, "3 

nuclear operator of the fluisi (see, for example, [12] ) t 


3c»j * 


C2.il) 






and o^ is the Pauli matrix; 


In these formulas " X ^ 

?!, and -- weak form fai^tors of the nucleon (for more details 
see the following section); p — impulse operator; M — nucleon 
mass. Summation is performed over all of the nuclei. Dots mean 


that the terms 




are omitted, and also corrections for the 


exchange meson currents. In addition, we have omitted the Induced 
pseudoscalar term (Included Ih the cross section only in terms which 
are proportional to the lepton mass) and it is assumed there are no 
currents of the second kind [l2l . 


We are only Interested in the cross section which is summed 
over all the spin and angular variables of the final nuclear system 
In this case, we may show (see Appendix A) that the convolution 
n n T. is reduced to the form: 


^ -y * » I 

Here 0 is the angle between t and ('scattering angle), * 

In deriving (2.5), we use a system of coordinates, in which the 


vector q is directed along the E-axls, so that q® q^,. The upper 
sign corresponds to a reaction with a neutrino ^ the lower — with 
an anti-neutrino. 

As was shovm in the Introduction, for calculations of the cross 
section we use a single particle shell model of the nucleus. The 
nucleus wave function is derived by the unique Slater determinant; 
the single particle wave functions of the initial and final states — 
the corresponding, functions of the Woods-Saxon potential belonging 
to the discrete and continuous spectrum. This model may be immedi- 
ately applied only to nuclei with filled J -shells. Therefore, in 
deriving the formulas for the cross section, we only consider these 
nuclei . 


In the single particle model, the components of the vector 
, as is readily found from C2.3), are calculated wxth the 

formula : 


^ { -iix A 

(<iV’ yi* c (x) t. C*) , 


( 2 . 6 ) 


where and are the single particle wave functions of the ini- 
tial and final states; — operator, 

+ ••• 

(2.7) 


The equation (2.7) includes only the terms corresponding to the 
components left in (2.^. In addition, we use the additional as- 
sumption v/c « i , where v is the velocity of the initial nucleon 
in the nucleus. In the nuclei ^/c <^0-25, so that the correspond- 
ing error is less than 25^. Due to this approximation, in (2.7) 
there are no terms with a derivative with respect to x, and the 


D • -F<i -Lf' 111 

’ iM ^ iM 


n . p - c ® ^ . 

*^•9 **• 






calculation includes radial integrals of only the simplest types. 
The validity of this approximation has been studied, for example, 
in [13], where it was shown that the discarded terms are only im- 
portant in a differential cross section for small scattering angles 
0 < 5°. 

A 

Following [1^], we shall use several modified operators . 

We shall replace the first component in (2.7) as follows; 

so that the vector part of the operator 3 satisfies the law of 
conservation of the vector current; 

Let us now calculate in our model of the cross section the 
reactions (1.1) summed over all the possible energy states of the 
final nuclear system. 


a) Transitions in the discrete spectrum 


In this case the integral (2.6) has the form: 


( 2 . 8 ) 


The convolution is calculated directly using formula 

(2.5). The exponent in (2.8), as is customary, is expanded in series 
in terms of spherical functions, making it possible to separate the 
radial and angular variables. As a result, we obtain the following 
expression (details of the calculation are given in Appendix B) : 
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Tfte notation in (2.9) has the following meaning: 

Q . S C^V\_(<^x)R. (jx,)y}dx 

where ri»^,e; is the single particle radial wave function^ 


( 2 . 10 ) 


— spherical Bessel function: 


ilk c L / 


( 2 . 11 ) 




Here are the Klebb-Jordan coefficients. Determination 

of the function ^ is given in Appendix B (except for T and C, 

they also depend on the quantum numbers £■ , ] of the initial and 
final states). 


Let us now use formula (2.1). We shall assume that the response 
Includes the nucleus as a whole (the effective mass of the target 
equals the nucleus). We thus find (we assume the leptons are 
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relativiEsfclc) ; 



/IS 


k, 


k,-E 

1 1 v,v I 

M, 2- 


= k-E 


(2.13) 


Here E is the nucleus excitation energy (difference between the 
energies of the final and initial level); — nucleus mass. Sum- 
mation in (2.12) is designated by and designates summation by 
pairs and 


Let us consider the simplest particular case, when the nucleus 
has filled 1-shells (doubly magic) and we disregard the spin-orbital 
connection. Under these conditions, the expression similar to (2.9) 
for transition from the shells 1^ in 1^ has the following simple 
form: 


2 ^ 

m. he, /* r ^ 

* ) 


( 2 . 1 ^ 1 ) 


S r-\ X® s 








(2.10a) 
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b) TranBltlor. in eonblnuous spoofcra 


For these transitions, the integral (2.C) may be v/ritten in 
the form: 

A 




( 2 . 15 ) 




Here Is the wave function of the continuous spectrum for a 

nucleon with the enorp;y p^/2m and the spin projection (vje use 
the representation of the channel spin Vpe.jW — the 

corresponding: radial function^ jfip — the elp.en function of the 
operator;^e^^ ,/i-i •(?)) . The calculation of 

is carried out slmllaWy to the case of transitions in 


r r 

a discrete spectrum: 


• \i I vd-VU^f • E bj-0 E [(-t) 

)l r I i r U V >• 


iS-j') 


M ‘•mW ”U 


0 




♦ P. F* i 1“' a. ) "iVj/, J 


( 2 . 16 ) 




As previously, the functions H and B are determined by the 

formulas (2.11), and the radial integrals now depend on the Impulse 



(a. 17) 








(2.16) Includeo Gummafcion ovor the projection of the spin of 
the channel ^nd over the direction of the nucleon fllr,ht. 

X > 

'fherofore, v/e must Integrate only with respect to the modulus p. 


According to 





(2.1), the cross section equals; 



( 2 . 18 ) 


Deserihing '•■he cross sections of th5.s form, wo assume that the 
flight diroction of the lopton and nucleon are independent, i.e., 
we disregard the limitations on the final state v/hich arise due to 
the law of conservation of momentum. Strictly speaking, we would 
have to expand the wave functions of nucleons in the nucleus with 
respect to plane waves and to introduce the law of conservation of 
momentum explicitly into the calculation, which would greatly com- 
plicate the discussion. Up to a certain degree, we can compensate 
for this defect, by introducing a term considering the response in- 
to the energy -function, and assuming that this response is per- 
ceived by a free particle of the mass M (effective mass of the 


target) : 




(2.19) 


Here M. and M- are the masses of the initial and final (produced 

11 X- 

after the flight of the nucleon) nuclei; E — excitation energy of 
the final nucleus . We obtain the following expression for the cross 
section : 


d.V 


G 


isuvl asf 


o'^o 




( 2 . 20 ) 
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and the connection of p with and the scattering angle is' found 
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from the condition that the argument of the 5 -function equal zero 
(2.19). 


If we disregard the spin-orbital connection in the wave func- 
tions of the continuous spectrum, then (2.15) is reduced to the 
expression 



( 2 . 21 ) 


and formula (2.16) is greatly simplified: 



2. Neutral current 


a) Elastic scattering of neutrinos by a nucleus 


This process, which occurs without excitation of the nucleus. 
Is only possible due to a neutral current. Within the framework 
of the single particle shell model (single-determinant wave func- 
tion of the nucleus) for the vector we have the following 

expression. Instead of (2.6), 




(2.23) 


where are the single particle wave functions, and summation is 
performed over all nucleons of the nuc!^eus . We shall assume that 
the nuclear operator of the current is a mixture of an 

isoscalar and the third component of the iso vector. For example. 


/16 
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for ) WG havo the following, instead of (2.i|): 

C51 . 1 F-X) <f; ‘ fA,) ?, 

.(ist*.2o 4S'C2-io|.)l 


( 2 . 24 ) 


and slmil^ly 
tain 

for 5(x) • 

, Which will ] 

operator A 

(compare with 

r 


K 

-4 

3*- 


/\ 

1 ' 



ZM 


( 2 . 2 >, 


Calculations, similar to thc>;;e carried out above for the case 
of a charged current, lead to the following formula for the cross 
section of the elastic scattering of the neutrino by the nucleus: 


^0 


aa (r.f i, !i*sw| V* f ^ » 

Mt *■ 

1 1,1 


( 2 . 26 ) 


(-4 




( 2 . 27 ) 

+ F..F. ( -I + ^ + • ]f^*i^^V 


Here we use the notation: 
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j 1 . 

<iv«i“'lji'^i><.UM^c'ol)^M„^<iou)luj).<;i'oio|c'o) • 


/ 1 7 




\-i\‘ ''l(^iv^')(.v<K^'‘^*K“'^*) 

JWjr 


(2.28) 




MiM^ 


Summation is performed in (2.27) over all paj.rs of nucleons in the 
nucleus (including i = ra) . The form factors and P^^ in (2.27) 

have the following meaning: if i is a proton, then C * F f P * i 

- 0 1 ' 1** I M .A ^ ’ 

if i is a neutron, then “m/V 


Pormula (2.26) generalizes the expressions in the literature 
[l6] to the region of average energies. Por small (E^ ~ 20 MeV) 

energies ("astrophyslcal" energy region), in (2.27) we may confine 
ourselves to terms with 1 • "3 * 0 . In this case 

1 (2.29) 

« 16 k.k: cs" |[ f; a » f( (i- k)] _ 
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It may be seen from (2.29) that the isosoalar part of the stream’ /1_8 
leads to the coherent cross section (a » 



In particular, ' » - . T'in tile Weinberg-Salam model [ 173 . 

b) Inelastic scattering (excitation of the nucleus) 

A separate calculation is not required here, since it is com- 
pletely similar to the case of the charged stream. The only change 
consists of the fact that now the cross section is divided into two 
parts — the proton (p-p transitions) and the neutron (n-n transi- 
tions). Each part, for a corresponding replacement of the form 
factors, may be calculated according to the formulas given above. 


3. Calculation results 

In deriving the formulas for the cross sections, it was assumed 
above that the nucleus-target has completely filled j -shells, and, 
consequently, there are unldetermlnant wave functions in our model. 
In order to approximately consider the real filling of the shells 
with nucleons, each single particle cross section describing the 
transition between two shells is multiplied by the coefficient [18] : 

. ■ 

Here N. is the number of occupied positions in the j^-shell, and 
is the number of free positions in the j^-shell. In this case 
of discrete-continuous transitions with j ^-shells, the correspond- 
ing single particle cross section is multiplied by 


/1 9 


All of the form factors of weak interaction for nucleons in 
the nuclei were selected just the same as for free nucleons. In 
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the case of a charged stream, using the hypothesis of conservation 
of the vector stream |l9j , we identified the weak vector form factors 
and Pg of the nucleon with the isovector electromagnetic Dirac 
and Pauli form factors : 


-fCvC'v') . 


C3.1) 


Thus F^(0l •! , 2MFj(0) . Here, = 1.79j 

= -I .91 — anomalous magnetic moments of the proton and neutron. 

Por P«(0), we use the value -1.2^1 I203 . The dependence of the 
a 2 

form factors on q may be described by the dipole formula 

( 1 * , i = V.A . (3.2) 

F.' •F.'ZMP., ■' 

Here = 8^0 MeV I2l3 . In the case of a neutral stream, we repre- 
sent the form factors in the form of the sum of the isoscalar and 
Isovector parts: 

^ I (3.3) 


Here the index 0 corresponds to the isoscalar part, and the index 
1 — to the Isovector part of the stream. In the Welnberg-Salam 
model 

F“io)‘-W’-e„ , F^\o) . ^ 

f/CoV»o, F;(.)4F,toK 

2 

Por the Weinberg angle 6 , we use sin 0 = 0.3* Then we have the 

w w 

following for the transitions between the proton states: 

F,p = , ' 


h'fW »F.f(0)*2.HFip(,0) = -0,i ♦ 
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For transitions between the neutron states 


For calculations of the Inelastic process cross sections (1.1), 

(1.2) j we have written a program making It possible to calculate 
the wave functions of the discrete and continuous spectra and the 
eigen energies of the connected states using the given nuclear 
potential. We can also calculate the cross sections of the processes 
examined using the formulas given in section 2. Appendix C gives 
greater details about the selection of the nuclear potential and the 
numerical calculations of the wave functions and the Interaction 
cross sections. In the calculations of the discrete-continuous 
transitions, we have assumed an effective mass of the target which 
equals the nucleon mass. 

For specific numerical calculations, we selected the following 
nuclei ^^0, ^'^01 ^^Fe, "^^Ga, ^^Br. The neutrino reactions with 

the nuclei are Important both for accelerating experiments (with the 

^^0, nuclei), and for experiments with solar neutrinos 

37 71 8l 

(^^Cl, "^Ga, "^Br) . The neutrino Interaction cross sections with 

37 71 8l 

detector nuclei (such as ^ Cl, '"^Ga, "^Br) must be known in a wide 

range of the energy changes (0-1 GeV) in order to evaluate the 
background of the cosmic ray neutrinos in the experiments [ 5 ) 


/21 


(the background arises from the interaction of the cosmic ray 
neutrinos with the detector nuclei). On the basis of the formulas 
obtained for the cross sections and the program, we can estimate 
the background for an arbitrary detector. 

In addition, for reactions of neutrinos with the nuclei exam- 
ined, other authors have made calculations, whose results can be 
compared with ours [1-8]. 

As a result of the calculations, we obtain differential and 
complete cross sections of inelastic reactions (1.1) and (1.2) 
occurring due to charged (n-p-transltlons ) and neutral streams 
(p-p and n-n-transitions ) . Figures 1-12 give the basic computa- 
tional results. 

Below, for brevity, we shall designate the discrete-discrete 
transitions as g-g, and discrete-continuous transitions as . 

Figures 1-6 show the cross section of the reactions 

(l.l) and (1.2) using the formulas given in the second paragraph, 
for the nuclei ^^0, ^^0, ^"^01, ^^Fe, '^^Ga, ®^Br. The solid lines 
show the behavior of the sum of the cross sections ^ 
and of the transition; the dashed lines show the individual 
cross sections of the transition. The numbers 1, 2, 3 

designate the transitions caused by a charged stream (n-p) and 
the transitions due to neutral streams (n-n) and p-p). The be- 
havior of the cross sections of the reaction (1.2) is similar to 
the behavior of the cross section for the reaction with charged 
streams (.1.1), since the difference in these cross sections is 
basically due to a difference in the form factors. The total cross 
section of the n-n and p-p transitions is approximately 30^ of the 
cross section of the reaction (1.1) almost in the entire region of 
the change in E , except for the region < ^0 MeV, which is close 
to the threshold of the reactions examined. 

Figure 7 gives the total cross section of fhe reac- /2 2 

tlon (1.1) for the nucleus ^^0 with closed n and p shells, along 
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with the results obtained by other authors. In the < 80 MeV 
energy region, which is close to the threshold, where the cross 
section of the transition is much less than the cross section 
of the transitions, the cross section we obt-ained is similar 
to the total cross section of the transitions between all indivi- 
dual levels examined in [33 , obtained by adjusting the experimental 
data. In the > 80 MeV region, the cross section of the 
transitions predominates; therefore, our cross section is greater 
than the total cross section obtained In [33- The shell calcula- 
tions performed using the approximation of completeness [8] in 
the < 200 MeV region, lead to a larger value of the cross sec- 
tion than with ours, and at E^ = 300 MeV it is approximately 2 
times smaller than ours. It is known that the study [83 could 
not explain the experimental behavior of the cross section 12 

obtained at Argonne [22] (the theoretical values were less than 

2 

the experimental values by a factor of approximately 2-3 at q 
0.05 GeV ). The results of our study show that this divergence 
may be reduced, if we consider the transitions more exactly in the 
continuous spectrum state. At = 300 MeV our cross section ap- 
proximates the cross section of the reaction (1.1) for 8 free 
neutrons 

Figures 1, 3, 7 use the crosses to designate the cross 

section for the reaction (1.1) obtained in the Permi-gas model in 
1 73, where the matrix elements were taken from experimental data 
with respect to the y-capture. For light nuclei ( C, 0), the 
method used in f7] is worse than for average and heavy nuclei, and 
the agreement with our shell calculations is not good enough, where- 
as for ^"^Cl and ^^Fe cross sections, which we obtained within the 
framework of the shell model and the Permi-gas model, coincide in 
the examined region of neutrino energy. The cross section formu- 
lated from the ^^0 on the basis of the Ferml-gas calculations from /23 
[I 23 behaved similarly to the cross section from [73 . 


ep, 


12 

for 


Figure 8 shows the total cross section 
[the reaction (1.1)] and the cross section obtained for the same 
reaction in [1] . The divergence between our cross section and the 
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cross section from [13 Is great for small E,,, in contrast to the 

16 ^ 

cross section of the reaction with 0, which is apparently re- 

12 

lated to the fact that the C is a light nucleus having a cluster 
structure to which our shell model is only slightly applicable . 


Figure 9 makes a comparison of the total cross section for 
one neutron for the reaction (1.1) for all of the nuclei examined, 
with the cross section of the reaction (1.1) for a free nucleon, 
given in [12] . It may be seen that the difference in the cross 
section 6(Ed/lJ is great at < 100 MeV. For large energies 

the differences become small, and at E^ ■•= 300 MeV the cross sec- 
tions approximate a single-nucleon cross section. 

Figure 10 makes a comparison of the cross sections els'/ Ail 
of the reaction (1.1) obtained in different models for the case 
^^0 at = 300 MeV. 


Figure 11 gives the reduced differential cross sections 

for all of the nuclei examined and makes a comparison 

with the cross section del Ail for a free nucleon. In connection 

with these figures, we should note the following. In the low angle 

region (for fixed energy E^, of the corresponding region of small 
2 ^ 

q ), the experiment gives much larger cross sections than theory. 

As was noted in the study of Cannata and Leonardl [133 , in this 
region it is necessary to consider relativistic terms in the ex- 
pansion of the nucleon stream, and also to consider the collective 
effects in the nucleus. 


Figure 12 shows the quasi-elastic peaks a«^/aadEj, for the 
reactions (1.1) in the nucleus ^*^0, calculated on the basis of our 
model and the Fermi-gas model [12] . The difference in the positions 
of the quasi-elastic peaks in these models is due to the fact that 
there are no thresholds in the Fermi-gas model for the reactions 
(1.1) and (1.2), and the nucleons do not occur in the potential hole. 
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^ . Conclusion 


We have described the processes of inelastic and elastic 
neutrino-nuclear reactions in a wide range E within the frame- 
work of a single-particle shell model of the nucleus with the 
Woods-Saxon potential . This description represents a step forward 
as compared with the simple model of the Perml-gas, since the use 
of the potential with a finite depth makes it possible to describe 
the state of the continuum and the related states using the solu- 
tion of the Schroedlnger equation with this potential. The parame- 
ters of the potential were determined so that the energies obtained 
of the connected states and the thresholds of the single-particle 
transitions were similar to the experimental values. We calculated 
the suras of all possible transitions in the connected excited states 
and in the state of the continuum for nuclei of 

'■^Ga, "^Br in the neutrino energy range of ® ^ ^ 300 MeV. An ex- 

amination was made of reactions taking place due to the charged 
streams and due to neutral streams of weak interaction. It is 
shown that the cross sections of Inelastic reactions occurring due 
to the neutral streams represent approximately 30^ of the value for 
the cross section of reactions with charged streams. 

The method described in this article is universal, and may be 
applied to any nuclei, and gives good accuracy for nuclei with 
closed shells. Only the parameters of the Wood-Saxon potential 
are necessary, which lead to a system of levels of the nucleus which 
is close to the experimental one. 

In calculations of transitions in a continuous spectrum, we 
disregarded the possibility of a nucleon sticking in the nucleus, 
and it was assumed that the potential is real. Por small q, the /25 
probability of a nucleon sticking in the nucleus is great. With 
an increase in q, the probability of a direct flight of a nucleon 
participating In the reaction Increases. Por a, more exact representa- 
tion of the state of the final nucleus, calculations are necessary 
which consider the imaginary part of the optical potential, and 
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also the ealculafclon of the evaporation proGesseo of nuGloono from 

the excited nuclei. In the high energy region, the results obtained 

using the shell model are similar to the results obtained using the 

Permi-gas model. The situation is more complex in the threshold 

energy range, since there it is necessary to consider the residual 

interaction of nucleons in the nucleus. Not considering this fact 

leads to the fact that the cross section we obtain exceeds the 

real one by approximately a factor of two for nuclei with unfilled 
1 ? 

shells of the C typo, as would be expected from comparing our 
cross sections with cross sections obtained using a careful adilust- 
mont of the experimental data tl]/[3l. 

It is also necessary to consider the residual Interaction in 
order to find a more exact distribution of the energy probabilities 
for the excitation of the final nucleus. 
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APPENDIX A 


Det us give the basic sfcagen^ln deriving the formula (2.0) 
for the convolution of First of all, by direct 

calculation using (2.2) we obtain the expression 


■3vw[no 1 


N 4 

Icfi f 

' Uc 

v.v;) 


. ” V.,lc 

. [IC'] 

1 - — 


k' \ 

k \c' 
*•0 

U. 

k.'j 


(A.l) 
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It is more convenient to use the fact that (see 1 231 ) the 
cross section in which we are interested is summed over all spins 
and integr»ated with respect to the flight angles of all particles 
(except for a lepton). The unique vector, which "exists" after 
the summation, is the vector q. Therefore, all three-dimensional 
vectors and tensors (not lepton) existing in (A.l) may be expressed 
(after summation) by the components of the vector q and the tensor 
and Consequently, we may write (the 

bar above indicates the summation): 


3.1* - , A> 






(A. 2) 


(A. 3) 




(A. 4) 
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It followG from (A. 2) - (A JO 
Dtitution directly before the 


that vre may make the fol*’owlnc oub- 
oumination, l.e., in (A.l); 




(the functions B, 


(A. 5) 

C and D are found from (A. 3) - (AJh) by removinc 


the bar above) . 


Vie selected the coordinate system so that the vector q is dis- 
tributed along the Z-axis . Except for *q, we have two lepton vectors 
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(A.S) 


The first and second terms in (A.l) may be transformed consider- 
ing (A. 5) as follows: 


(^n- ■♦w ) Jo ^ ^ 2. 


(A. 7) 


The first component in the last term in (A.l) may be discarded, 
since ' r7 rl'TU-* 


2.JJ 




m ^ t \ 

[UU ] 


-rrr-r*o 


k.W ■ ' k.W<! 

/ 


(we recall that - k*'k. ), 

Finally, we transform the last component in (A.l): 
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% 
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ColleoMnB all the formulas, we reaeh the expreoslon 



(2.5) 


APPENDIX B 


We shall explain in general features how we derived the 
)rmulas (2.9) > (2.16) and (2.27). ¥e first perforin summation 
over the spins and perform integration with respect to the angles 
in the integrals (2.6). Expanding the exponent in this integral 
in series in terms of spherical functions, 
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(B.l) 


Vie reduce this problem to calculating the following two standard 
Integrals : 


if 1 1 I - jj- • < t.o 1 0 le,o) . 
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' jj 3 1, 

We also use the following formulas for ®j and ®j symbols [24] 
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function KiO is determined by the formulas: 

, , , - se, 

. 

Substituting (B.3) into (B.2), we obtain simple expressions for 


(B.4) 


Integrals : 
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The Klebb-Jordan coefficients, containing and are lost 

during squaring and summation over and M^, 


I s:„. . 


For example, we obtain the following result for the term in 
(2.5) j which is proportional to 



[O0«lco><*«|u>^^ 1 1 leo>f itK) |(Vc)| j 


The Interferential terms, which are proportional to the product 
of the vector and axially-vector form factors, are reduced during 
summation over T and T' always, with the exception of the last com- 
ponent in (2.5). 


We may readily obtain the formulas for the particular cases 
l(2.l4) and (2.22)] not by using (B.3), but directly by using the 
following well-known relationships in (B.2) 1253: 




k «■ ii 
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APPENDIX C 
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Selection of Effective Nuclear Potential and 
Certain Details of Numerical Calculation 


1. Discrete spectra 


The wave functions and the eigen energies of the discrete 
spectrum were calculated by solving the Schroedinger equation with 
the effective nuclear potential in the Woods-Saxon form (consider- 
ing the Coulomb potential for the protons). The potential was taken 
in the following form: 


V(:^^ V,Cx) -V.f.Cx) i . 

r \-i 


(c.l) 


t » 0 s • 





(C.2) 


The Coulomb potential Vq(X) was assimed to equal the potential of 
a uniformly charged sphere with the radius and the charge 


V.Cx) .(a-4) 1X^(3^), 


(C.3) 


where j for the protons 



for the neutrons 


at XiR. 
at X'vR. 


‘>5 ■ (C-5) 

To solve the Schroedinger equation, a program was written making it 
possible to find the wave functions and the eigen energies of all 
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the states of the discrete spectrum using the given A and Z. 
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10 1 

For the and 0 nuclei, we made several detailed calcula- 
tions with different values of the parameters for the nuclear poten- 
tial in order to obtain partial transitions in the reactions (1.1) 
and (1.2), which were similar to the experimental values for the 
thresholds. This was due to the fact that, as the calculations 
showed, the cross sections of the partial transition between two 
levels depend mainly on the difference between the eigen energies of 

thesb levels, i.e., on the threshold, and to a lesser extent, the 

12 

asbsolute values of the energies of the levels. For C, we final- 
ly selected the following set of the potential parameters (0.1) - 


(0.5) : 


Mev V 


? 




Vv 




Mev 




( 0 . 6 ) 


The quantities 
the formula 


fto and R.,, 


were calculated in the case 


12 


C using 




(0.7) 

For the nuclei ^^0, ^'^01, ^^Fe^i "^^Ga, ^^Br, the values of ^ 
and Vt for the potentials (C.l) - (0.2) were taken on the basis of 
the study [261. These parameters are given in Table 1. 


Ing: 


The values of the quantities V ^ and V for ^^0 are the follow- 

^ on op 


V 


Otw 


= 56 MeV Vop = 52 MeV , 


(0.8) 


Q Y 71 8l 

For ^'01, ^°Fe, "^Ga, "^Br, we used the interpolation formulas from 

[27] for and 


Vo^= 49,84 46,5* (MeV) 

57,8 + 9,2 • (t'i-'Z.VA 
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The quantities R* and Rj foi-* these five nuclei were determined 
using the formula [26] : 

, t«o,s. (c.io) 

The potential (C.l) - CC.5) with parameters determined from (0.9) 
and (0.10) and Tfacle 1 make it possible to obtain the energy of 
the upper filled single-particle levels, coinciding with the ex- 
perimental data in a wide range of A and Z. Table 2 gives the 
values of single-particle energy levels for all of the nuclei which 
we examined, obtained as the result of numerical calculations. 


The Schroedlnger equation for the radial part of the wave func- 
tion the nucleon, located in a field with the potential 

(0.1) - (0.5), was solved in the usual way [28] using the substitu- 
tion of C*^) 0^) j as the result of which the following 

equation was obtained for 'i -vf, W 

/.ijW - W .v.f» V4 i 


ICx) 




( 0 . 11 ) 


where ^0 is the eigen energy of the level with the quantum 


numbers W;<,L j j . 

The solution (O.l) must satisfy the boundary conditions 
/ V /• \ /\ and the normalization condition ( x J 


In a small region close to the center of a nucleus 
0 <C X ^ ^ the solution is selected in the form of the spherical 

Bessel function 




( 0 . 12 ) 
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where for the 
tons 


neutrons 


Xa was selected to equal 0.4 f. 


and for the pro- 
The quantity 


As the solution at X oo , we selected the Whittaker function 

. "here K. ^ {a ^ y 0,0 

for protons and 0 for neutrons. Using its Integral repre- 

sentation, we may write 





r 


^ dt 


(C.13) 


In the case of the neutrons *^ = 0 and the function W is re- 
duced to the Hankel function of the imaginary argument, which can 
be represented In the form of the finite series: 


where O » K.X • 
factors. 



The quantities N 




(O.ll)) 

Ng,N^ are arbitrary numerical 


In the case of protons, we obtain the function W and its deriva- 
tive with respect to x, calculating the Integral from. (C.13) and the 
integral of a similar type arising in the calculation of the deriva- 
tive , using the Laguerre formula with 6 nodes [291. 

Equation (C.ll) was solved using a computer as follows. We 
assigned the quantum numbers 'and the value of ^ j 

From. 0 to the point x^, the analytical expression (C.12) was used 
as the solution. The values of the function (C.12) and its deriva- 
tive with respect to x of the first order were calculated at the 
point x^j and were used for a numerical solution of the complete 
equation (C.ll) using the Runge-Kutta method of fourth order with 
an automatic selection of the Integration steps, with the condition 
that at each step the relative error in the value of the wave func- 




tion and Its derivative did not exceed 10”^. The numerical integra- 
tion was performed from the point x^ to a distance of 13 fm from the 
nucleus center. By roughly varying the energy , it was found 

that the number of nodes of the calculated wave function became equal 
to n+lj including the node at the point x = 0. Then, using the cal- 
culation of the relative deviation of the derivative of the wave func- 
tion, obtained by numerical integration from the derivative of the 
asymptotic Cc.l3) or CC.l4), we found a more exact value of the 
eigen energy. The calculation of the eigen energy of the level was 
completed, if the relative deviation of the derivatives at the point 
where it interlaced with the asymptote, which was a distance of 
13 fm from the center of the nucleus, was less than 10“^, or if the 
absolute error in the value of the eigen energy was less than 
5 '10 MeV. Then the wave function obtained was normalized to 


unity. To find one value of the eigen energy, from 15 to 23 itera- 
tions were necessary. The values of the wave function corresponding 


to the values of the argument 13 x^, where x^ were the nodes of the 
quadratic formula ^ ^ _ j_ ^ 1 . M 


were calculated and written on magnetic tape, since they were neces- 
sary to calculate the integrals 0 . . and (? 


from expressions (2. IQ) and (2.17). 
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To correct the program, we calculated the wave function and j 

the eigen energies of the nucleus of 208p^ with the potential 'i 

parameters from the study of Blomqvlst and Wahlborn 128]. The re- i 

suits obtained fully coincided with, their results. j 

1 

2. Continuous spectrum i 

The states of the continuum were described by wave functions ^ 

obtained by solving the Schroedlnger equation with the same poten- -j 

tial which was used to calculate the wave functions describing the | 

connected states. One difference was that, to reduce the compute- 1 

tional time, the spin-orbital term in the potential was discarded. I 

'I 

In addition, we only calculated the wave functions with t 6 10. i 
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The correction showed that the states with large ^ did not make 
a great contribution to the cross section in the examined region 
of the neutrino energy. 

The equation for radial parts of the wave functions of the 
continuous spectrum of" ^ nucleon with the energy Ej 

after substitution of W ^ ^ was obtained from 

(C.ll), if we make the substitution “'E and discard the 

J 

spin-orbital term in the potential: 

-O.ObWti-O'^Sfi.W • 

In the asymptotic region from (C.5) we obtain the equation for 
Coulomb wave functions FrW and [15]. 

Close to the center of the nucleuSj (C.15) has the solution (C.12) 
where the expression for 0 is obtained from the previous substitution 

**E. 


The expression (C.12) was selected as the wave function in the 
region from 0 to the point x^^, and from the point x^ to a distance 
of 13 fm from the center of the nucleus the equation (C.15) was 
solved numerically by the Runge-Kutta method with a constant step. 

The step was chosen so that their lengths did not exceed 0.1 fm, 
and was constant In the Intervals of ^ where x^ are 

the nodes of the quadratic formula, ]ust as in the case of the wave 
functions of the discrete spectrum. There had to be an integer 
number of integration steps in these intervals. When these condi- 
tions were observed, the Integration accuracy was sufficient, and 
the values of the function calculated at the points 13 x^, after 
they were multiplied by the number m, determined from the interlacing 
condition with the asymptote, were written on magnetic tape. A dis- 
tance of 13 fm from the center of the nucleus, the numerical solu- 
tion was Interlaced with the asymptote, which was selected in the 
form [ 15 ]: 
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■ V Iff [u%s;-Pt w « ■ G-eCx)i. 


(c.l6) 


Since we have not considered the imaginary part of the potential, 
the phase of the scattering S*| is real. The quantity m , which 
must be multiplied by the wave function, obtained by numerical in- 
tegration, and the phase are determined from the Interlacing 
conditions : 

( 0 . 17 ) 


VW * . (0.18) 

3 . Calculation of integrals and special functions 

The Integrals ^ and j. 

included in expressions (2.10) and (2.17) were calculated. The 
Integral from 0 to 00 was represented in the form of the sura of 
two integrals, and the first had integration limits from 0 to 13 fm. 
It was calculated using the quadratic Gauss formula with l4 nodes, 
using the values of the wave functions written on magnetic tape. 

The second Integral with limits from 13 fm to 00 was calculated 
using the asymptotic expressions (0.13), (C.l4) and (C.I 6 ) for the 
wave functions. Since expression (0.13) contains the exponent with 
negative index, this integral may be readily transformed to the form 
, after which it is calculated using the quadratic 
Laguerre formula with 5 nodes. 

The Coulomb wave functions , the spherical Bessel functions , the 
Klebb-Jordan coefficients used in the study, were calculated using 
specially written programs, whose correctness was verified using 


well-known tables of corresponding values. 
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HEADINaS FOR THE PIOUREB 


/A 2 


Pig. 1 . 


Pig. 2 
P16. 3 
Pig. ^ 
Pig. 5 
Pig. 6 
Pig. 7 


Pig. 8. 


Pig. 9. 


Pig. 10. 


Total cross section of reactions (1.1) and (1.2) for 
The solid lines show the sum 
of total cross sections and tJH positions. The 

numbers 1,2,3 designate the cross sections of the n~p, 
n-n and p-p transitions. The dashed lines show the 
cross sections of the gg transitions. The crosses show 
the cross section for the reaction (1.1) according to 


17). 


16 , 


The same as for Pig. 1, but for the nucleus 0. 


The same as for Pig 
The same as 
The same as for Pig 


1, but for ^'^Cl. 

for Pig. 1, but for ^^Pe . 

71 

1, but for ' Ga. 

81, 


The same as for Pig. 1, but for Br . 

Comparison of the total cross section 

16, 




of the 


reaction (1.1) for *’^0, obtained in this study, with the 
results of other authors . The numbers designate: 

1 - the study i 2 - cross section from [ 3] i 3 - cross 
section obtained in the Permi-gas model [12] ; 

4 - cross section from £sj ; 5 - cross section for a 
free neutron [12] , multiplied by 8^ the crosses — cross 
section according to [7] . 


Comparison of the cross section 


O'- 


(solid line) 




and ^ (dashed line) of this article with the transi- 
tion cross section considering all individual levels from 
[11 (dot-dashed line) . 

Comparison of the reduced cross section 
with the cross section for a free neutron [for the reac- 
tion (1.1)]. . 

Comparison of the angular dependences dS’/JllL for ^^0 in 
in different models: 1 - this article; 2 - Permi-gas 

model Il2]; 3 - 18]; 4 - free neutron. The graphs are 
compiled for = 300 UeY [the reaction (1.1)1. 


/43 


39 


Pig. 11. 


Reduced differential cross sections 
for all nuclei examined and d^liSL 
= 300 Mev. 


for a free neutron j 


Pig. 12. Quasi-elastic peaks cI^MIIoIE for the reaction with 

1 6 ^ 

0. The solid lines show the quasi-elastic peaks ob- 
tained in this article; dashed lines — In the Permi- 
gas model [12], The numers around the curves designate 
the quantity , Where 9 is the angle between 

the neutrino and the electron. The neutrino energy 
equals 300 MeV. 
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